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ABSTRACT 

We examine the distribution of transition discs as a function of mm flux. We confirm 
that as expected in any model in which most primordial discs turn into transition discs 
and in which mm flux declines with time, transition discs have lower mm fluxes on 
average than primordial discs. However, we find that the incidence of transition discs 
does not, as expected, fall monotonically towards large mm fluxes and we investigate 
the hypothesis that these mm bright transition discs may have a distinct physical 
origin. We find that mm bright transition discs occupy a separate region of parameter 
space. Transition discs in the bright mm sub-sample have systematically higher ac- 
cretion rates than those in the faint mm sub-sample, along with being systematically 
weighted to earlier spectral types. 

Key words: planetary systems: protoplanetary discs - stars: pre-main-sequence. 



1 INTRODUCTION 

Protoplanetary discs with apparent opacity deficits at Near- 
IR (NI R) waveleng t hs we re fi rst identified in early sam- 
ples bv lStrom et~afl <|l989h and lSkrutskie et al.l (|l990h . The 
lack of protoplanetary discs observed to be glob ally opti- 
cally thjri_at_ajf_j^ : dii_(e ; g. [e l201ll : but see 
also ICurrie fc Sicilia-Aguilarl 120111 ). led to the interpreta- 
tion that, if inner-hole objects indeed represent an evolu- 
tionary stage from a primordial, optically thick state to a 
disc-less state, then disc dispersal must preferentially oc- 
cur from the inside-out. Under this assumption they were 
consequently christened 'transition discs', and their rar- 
ity allowed the calculation of a disc clearing time-scale of 
around 10 — 20% of a disc's primordial lifetime for both solar 
type stars (e.g. iKenvon fc Hartmannl 1 19951: iLuhman et all 
l20ld ), and low mass stars (e.g. lErcolano et alj|2009t l201ll) . 
These transition disc observations have led to numerous 
models which attempt to explain their origin. One cate- 
gory of models explains this opacity deficit in terms of 
clearing of both gas and dust, for exam ple by: planet 



formation ( e.g. lArmitage fc Hansenl Il999l ); photoevapo- 
ration (e.g. IClarke et al.l 120011 ) and MRI driven winds 
(|Suzuki fc Inutsukal2009l ). The ofther category of models in- 
stead posits only a deficit of dust o pacity in the inner disc, 
resul ting from: grain growth (e.g. iDullemond fc Dominikl 
120051 ): dust filt ration by planets (e.g. Rice et alT 20061 ) or 
photophoresis jKrauss et al.l [2007) . At this stage it is not 
clear observationally or theoretically whether one of these 
processes dominate, whether they are all occurring, or tran- 
sition discs are caused by some yet unidentified mechanism. 

The observed sample of transition discs has now grown 
to a sufficiently large size that one can b egin to comp a re the 
prop erties of these discs to the models. I Owen et al.1 Hon], 
|2012| ) compared the predictions of X-ray photoevaporation 
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to the sample of observed transition discs in the -Rhoie and 
M* plane and identified that only approximately ~ 50% 
were consistent with being created through photoevapora- 
tion: namely, those with small holes iihoie < 20AU and low 
accretion rates M* < 10 -8 Mq yr _1 . Furthermore, several 
other authors have used properties such as the surface den- 
sity profile, presence of gas and/or dust in the inner hole, or 
disc mass to postulate the o rigin of individual transition 



discs (e.g. Cieza et all 20081 . l2010l : lAlexander fc Armitagel 



120091 : lEspaillat et al.|[2O10l ) 



Perhaps the most surpris i ng rec ent development is that 
presented by lAndrews et al.l (|201ll ): looking at discs with 
large holes (and therefore imageable in the sub-mm) they 
calculated that ~ 25% of discs in the upper quartile of the 
sub-millimetre (mm) flux distribution for Class II discs may 
be identified as transition discs. This result is contrary to 
our inference that transition discs represent an evolutionary 
stage where discs are transitioning from a primordial to disc- 
less state at the end of th eir lifetimes, as mm flux is observe d 
to decrease with age (e.g. lAndrews fc Williamdl2005l . l2007l ). 
In this scenario one would have naively expected a transition 
disc fraction that decreases monotonically with increasing 
mm flux. 

The analysis of statistically significant studies of ob- 
served transition disc samples has only recently become pos- 
sible, allowing us to investigate their properties and iden- 
tify possible sub-samples, taking care not to throw away the 
valuable information contained within upper-l imits (using 
Survival statistics e.g. lFeigelson fc Nelson|[l985l ). In this let- 
ter we use a large sample of transition discs (not just those 
which are imageable in the sub-mm) collected from the liter- 
ature, to investigate the ratio of transition discs to primor- 
dial discs as a function of the full mm flux d istribution. We 
can t hen compare this sample to both the (|Andrews et al.l 
l201lf ) result and our naive theoretical assumptions as to 
what it should be (i.e. peaked a low mm flux and then mono- 
tonically decreasing). In Section 2 we describe the transition 
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disc sample and look at its basic properties. In Section 3 
we discuss the more detailed properties and possible sub- 
samples of the transition disc population and in Section 4 
we draw our conclusions. 

2 SAMPLE 

We have collected a large sample of transition discs from 
the literature that have measured mm fluxes; we choose 
to use 1.3mm fluxes instead of 880/im fluxes as the for- 
mer are more readily available for the collected data. In the 
case that only 880/^m or 1.2m m measurements are available 
we use equations presented in ICieza et al.l (120081) fused to 
calculate disc mass - based on lAndrews fc Williams! [20051 . 
120071 ) to convert measurements into 1.3mm fluxes, and then 
we normalise all fluxes to a common distance of 140 pc 
using the Spitzer c2d distanc es. The transiti o n dis c s use d 
in the sampl e are t a ken from: Andrews et al.l (|201ll, 2012h; 



(200d 


2010h; Espaillat et al. 


J2001T 


20101); [Hughes et al. 


(201Q 


); iKim et alj ( 20091) ; iMerin et al. 


(2010h;lNaiita et al. 


(2007^ 





Many of these articles discuss the same transition discs, 
so we take care to remove multiple instances of the same 
source. In addition we remove any 'transition' disc which is 
caused by a known close stellar companion (e.g. CoKu-Tau 
4 Ireland & Kraus 2008). The total sample includes a to- 
tal of 76 objects classified as transition discs by the various 
authors. For the 76 sources that make up our sample, 60 
sources have available accretion rate and spectral type de- 
terminations; additionally 32 sources have had inner holes 
measured either through spectral type fitting and/or mm 
imaging. 

It is worth noting that lCieza et all (|2008l . I2010T I do not 
fit their sources for inner holes, and they use a transition 
disc definition th at is based on mea surements of the Spitzer 
colours. However, ICieza et all (|2010h do list a 'turn-off wave- 
length' representative of the dust temperature for the in- 
ner hole; thus we follow lOwen et al l (|201l| , |2012l ) in turning 
these wavelengths into conservative upper limits. 'Turn-off 
wavelengths' of 8/im and < 5.8/im are assigned inner hole 
radii of < 10AU and < 5AU, respectively. This increases the 
number of objects with any kind of inner hole measurement 
to 50 objects. We ca ution in the case o f the four sources 
that o verlap between ICieza et al.l ([2010h and iMerfn et al.l 
l[2010! ) , only two of these sources are identifie d as transition 
discs with inner holes by IMerfn et all ([2010T I. Therefore, it 
is importan t to bear in mind that the sources classified by 
ICieza et al.l ([20081 . |2010| ) may not all be discs with inner 
holes, and the only definitive method of inner hole measure- 
ment is imaging. 

2.1 The millimetre distribution 

We use the transition disc sample to construct the Kaplan- 
Meier product-limit (KM) estima tor (an estimate of the sur - 
vival probability P(xt > X) e.g. iFeigelson fc Nelson! 0*9851 ') 
for the mm distribution of transition discs. Furthermore, to 
compare the transition disc sample t o a sample of primor- 
dial discs we use the sample of Taurus ( Andrews fc Williams! 



l2005h; Ophiuchus (lAndrews fc Williams! 120071); Chameleon 



( Henning et"ai1ll993l ) fc Lupus ([Nuernberger et al.lll997T l to 

construct the KM estimator for the mm distribution of class 
II (transition discs have been removed) objects. Since our 
class II sample is exclusively from young (l-2Myr) clus- 
ters we further restrict our comparison to transition discs 
in the following young clusters: Taurus, Ophiuchus, Ser- 
pens, Chamelon & Lupus, as the mm dist ribution of class 
II discs does appear to evolve with time (|Lee et al.l l201ll : 




Figure 1. The survival distribution for the mm fluxes of the pri- 
mordial discs constructed from the young class II objects, — thick 
black line. Young transition discs are shown as the thick dashed 
blue line, while the dotted line shows the same transition discs 
sample excluding the Cieza et al. objects, and the thin dashed 
line shows the total transition disc sample. The error bars show 
the la errors in the KM estimator. The thin lines represent the 
smoothed distribution functions described in the text. 
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Figure 2. Figure showing the relative ratio of transition discs 
to primordial discs in each quartile of the underlying primordial 
mm distribution, constructed from the KM estimator (solid) and 
smoothed distribution functions (dashed). The values are scaled 
such that the fraction in the lower quartile is unity for the KM 
estimator. Error-bars show representative errors estimated from 
the lfj errors in the KM/KS estimators. 



iMathews et alj|2012l ). This leaves us with 48 transition discs 
to compare with our class II sample, although we note the 
mm distribution of these 48 discs is statistically indistin- 
guishable from the total sample (dashed line in Figure 1), 
thus we use the total sample later when looking at the prop- 
erties of transition discs. 

We show the two survival distribution functions for 
the transition discs (thick dashed line) and primordial discs 
(thick solid line) with l<r errors in the KM estimator shown 
as error bars in Figure [T] We note that removing the Cieza 
et al. transition discs - all of which may not be discs with 
cleared inner holes, as their SED's were not fit for inner holes 
as discussed previously - from the sample (dotted blue line in 
Figure 1) does not significantly affect the resulting survival 
distribution. As survival statistics does not account for the 
fact astronomical meas urements contain errors i n measure- 
ments and upper limits ([Feigelson fc Nelson!! 19851) which can 
mask an observational cut-off in a luminosity functions with 
true cut-off in a luminosity function, we compute smoothed 
distribution functions by broadening each measurement or 
upper limit with an appropriate Gaussian or error function, 
allowing us to make sure a hard observational cut-off is not 
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driving our conclusions seen in the data These smoothed 
distribution functions agree with those constructed from the 
KM estimator, indicating the use of exact upper limits in the 
survival analysis is not significantly affecting the result and 
our results are representative of the true luminosity func- 
tions. 

A glance at Figure[T]clearly shows that the mm distribu- 
tions of primordial (class II) and transition discs are differ- 
ent, and the hypothesis that they come from the same under- 
lying distribution can be quantitatively rejected at a level of 
5.7cr and 6.4(7 using the Gehan and logrank hypothesis tests, 
respectively. Furthermore, we can use the smoothed distri- 
bution functions to perform a Kolmogorov-Smirnov (KS) 
between the two populations, where we can reject the null 
hypothesis at the level of 2.4a, again indicating the use of 
survival statistics is appropriate for this data set. Addition- 
ally, we can use these two distributions to construct the rel- 
ative fraction of transition discs to primordial discs in each 
quartile of the class II mm distribution function, using both 
the KM estimator and the smoothed distribution function. 
This result is shown in Figure [2] where the relative fraction 
is scaled to unity in the lower quartile. 

Figure [2] shows a clear excess of transition discs in the 
lower quartile of the mm distribution, and a distribution 
function that is inconsistent with one that monotonically 
falls with increasing mm flux. Note that the transition disc 
and primordial samples do not come from the same star 
forming regions. If we sub-sampled the transition discs such 
that they are also in the come from the regions with com- 
plete class II mm samples, the sample size would be too low 
to be statistically useful. Therefore, Figure [2] only contains 
information on the relative ratios of transition discs to pri- 
mordial discs in different quartiles of the primordial disc mm 
flux distribution. 

3 DISCUSSION 

The high fraction of transition discs at large mm fluxes is dif- 
ficult to explain with the canonical view of transition objects 
as discs caught in a final evolutionary stage before dispersal. 
In this picture one would expect a steeply decreasing tran- 
sition disc fraction with increasing mm flux, as the latter 
also decreases with time. The high transition disc fraction 
at large mm flux contradicts this simple view and can be 
interpreted as evidence that either: (i) 'transition discs' are 
a misnomer and these objects do not physically represent an 
evolutionary stage from disc-bearing and disc-less systems; 
or (ii) discs with inner holes can have multiple physical ori- 
gins and the two populations shown in the data indicate 
two distinct classes of objects. The latter hypothesis will be 
further investigated in the following sections. 

3.1 Spectral type distribution 

Perhaps the simplest interpretation of the observed distribu- 
tion of transition discs with mm-flux is the existence of two 
population of transition discs, one peaked at low mm fluxes 
which would be consistent with the classical picture of tran- 
sition discs i.e. disc clearing, and a population of inner-hole 
objects with large mm fluxes which has a different physi- 
cal origin. If we split the transition disc sample into two 
sub-samples about the median Class II mm flux, 42 and 34 
objects fall in the low- and high-mm flux sub-samples, re- 
spectively. Discs with upper limits greater than the median 
value are placed in the sub-sample of discs with high mm 
fluxes, although we note this only represents a very small 
fraction (4 objects) of the total sub-sample. 

We can test the hypothesis that the mm faint transi- 
tion discs represent a later evolutionary phase in the life 
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Figure 3. Figure showing the distributions of spectral types for 
the primordial discs (solid), mm faint transition discs (dashed) 
and mm bright transition discs (dotted). The error bars indicate 
lcr errors. 



of primordial discs by checking that they share the same 
distribution of spectral types as the primordial disc popula- 
tiorE The similarity in spectral type distribution between 
the mm faint transition discs and the primordial disc sam- 
ple is demonstrated by the dashed and solid lines in Figure 
3: a KS test reveals that the two distributions have a high 
probability of being drawn from the same population (we 
can only reject the hypothesis they are different at ~ 0.6<r). 
On the other hand, the distribution of spectral types for the 
mm bright sub-sample of transition discs is clearly different 
from that of the primordial population and the hypothesis 
that these are drawn from the same population can be re- 
jected at the 4.8(7 confidence level. Clearly the bright mm 
sample contains a much higher proportion of stars of earlier 
spectral type than the primordial (or mm faint transition 
disc) sample. 

This finding needs to be interpreted with care. The 
primordial sample contains a large sample of objects from 
Taurus-Auriga (~ 40%) which is well known (e.g. Goodwin 
et al. 2004) to be deficient in stars of earlier spectral type. 
On the other hand, the transition discs are a serendipitous 
sample culled from clouds at a range of distances and there 
can be expected to be an obvious over-representation of sys- 
tems of earlier spectral type since these are more luminous. 
However, it should be stressed that these biases do not affect 
the obvious difference in spectral type distributions between 
the mm bright & faint transition discs since mm flux forms 
no part of the object selection and identification. 

It is also worth noting that although the mm bright 
transition discs contain a high fraction of earlier type stars, 
there is also a significant population of mm bright discs at 
late spectral types (see Figure 5). In fact, if we construct 
a figure similar to Figure 2 but include only stars of spec- 
tral type later than both T c g = 5000 & 4500K we obtain a 
qualitatively similar result, with a concentration of sources 
in the lowest quartile and a flat or gently rising distribution 
in the other quartiles. 

3.2 Transition disc properties 

We now examine the derived properties of transition disc in 
the two sub-samples. Figure 2] shows the distributions of ac- 
cretion rates. It is clear that distributions of accretion rates 
are different for the two transition disc sub-samples, with 

1 This test also assumes that there is no significant variation in 
transition disc lifetime across the limited range of spectral types 
(~M to ~G ) in the samples. 
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Figure 4. The distribution function of accretion rates for mm 
bright transition discs (dashed) and faint mm fluxes (solid), the 
nominal accretion rate detection limit of 10 -11 Mq yr _1 is shown 
as the dotted line. 



the mm bright transition discs having considerable higher 
accretion rates. We confirm the difference in the distribu- 
tions at the level of ~ 4.6a and ~ 5<r for the accretion rates 
using the Gehan and logrank hypothesis tests, respectively. 
We do note that the mm bright sample does appear to have 
in general larger inner hole radii (see Figure 5) , although as 
the majority of inner holes are not measured or fit for the 
mm faint sample we do not draw any statistical conclusions. 

Along with looking at the individual distribution func- 
tions for inner hole statistics we can look for correlations be- 
tween the transition disc statistics and the parameter space 
they occupy in the individual sub-samples. It is of particular 
interest in order to test the models to study the distribution 
of observed transition discs in the 7?hoie--M* plane as shown 
in the upper panel of Figure [5] We see that, broadly speak- 
ing, the mm faint objects (filled symbols) occupy t he region 
of parameter space identified bv lOwen et ail (|2012l ) as being 
conceivably produced by photoevaporation whereas almost 
all the bright mm objects (open symbols) lie outside this 
region. We also (lower panel) consider the relationship be- 
tween inner hole radius and stellar effective temperature in 
the two sub-samples. We note that in both panels it is im- 
possible to draw any conclusions about correlations between 
these properties within the faint mm sample on account of 
the large number of upper limits in hole size and accretion 
rate: a large f raction of our faint mm sample is drawn from 
the sample of lCieza et al.l (|201Cll ) for which we have crudely 
estimated limits on the hole size (see Section 2). 

However, correlations do appear in the mm bright tran- 
sition discs, where there is clear evidence of a correlation be- 
tween inner hole radius and accretion rate (we can reject the 
null hypothesis of no correlation at the level of 2.4<r), with 
large holes generally havi ng higher a c cretion rates. Such a 
correlation was noted by iKim et al.l |2009t) with a much 
smaller sample o f transition discs, and it is confirmed here. 
IKim et all |2009l ) also noted correlations between accretion 
rate & inner hole size with properties of the underlying star. 
The lower panel of Figure 5 shows clear evidence for a corre- 
lation between inner hole radius and spectral type (a crude 
proxy for stellar mass) in the mm bright sub-sample (we 
can reject the null hypothesis of no correlation at the level 
of 2.5<t). Moreover, it also clearly shows that the holes in 
the mm bright discs are systematically larger than those in 
the mm faint discs at all spectral types and therefore the dif- 
ference in inner hole radii shown in the left hand panel of 
Figure 4 is not purely a result of the mm bright discs being 
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Figure 5. Figure showing correlations between transition disc 
properties. The upper panel shows inner hole size and accretion 
rate, the lower panel shows inner holes size and effective stellar 
temperature. Triangles indicate the upper limits in inner hole ra- 
dius, squares indicate actual measurements. Open symbols arc for 
mm bright transition discs, and filled symbols are for mm faint 
transition discs. Objects at accretion rates of 10 -11 Mq yr _1 arc 
upper limits in accretion rates. We note some symbols contain 
multiple objects. The shaded r egion represents th at predicted by 
photoevaporation; tak en from lOwen et al.l l|2012h - note the ver- 
sion of this Figure in lOwen et al. I | |2012|) contains more objects 
that do not have available mm fluxes. Points with directly im- 
aged holes are also circled. 



of earlier spectral type. We note in passing that the correla- 
tions are still present (non-zero correlation coefficient) when 
one only considered those discs where the inner hole radius 
has been measured by imaging although it is less statistically 
significant where we can only reject the null hypothesis of 
no correlations at the 1.3er level in both cases. Finally we 
note that there is no evidence for correlations with mm flux 
for the inner hole radius, accretion rate and spectral type in 
the two individual sub-samples. 

We speculate the interpretation of these results is that 
transition discs are not dominated by one physical mecha- 
nism, but there are perhaps two mechanisms at play. One 
producing transition discs in mm faint discs (likely associ- 
ated with the end of the discs life) with low accretion rates 
and small holes, with the most likely origin being X-ray pho- 
toevaporation (Owen et al. 2011). The second mechanism 
may not be associated with the end of the disc's life at all, 
or the processes itself cause a high mm flux; the high accre- 
tion rates and presence of small dust grains inside the mm 
cavities (Dong et al. 2012) point to planet formation (e.g. 
Rice et al. 2006) . However, we caution that much theoretical 
work must be done to reconcile these observations with the 
planet formation scenario since - in addition to reproducing 
the correlations suggested in this work - there is consider- 
able fine tuning required in the dust distributions that can 
match both the SED and scattered light/mm images. 
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4 SUMMARY & CONCLUSIONS 

We have used a sample of transition discs collected from the 
literature to look at their statistical properties, and compare 
their mm properties to those of primordial discs (Class II 
objects) in Taurus Auriga. We demonstrate that transition 
discs are concentrated towards the lower end of the mm flux 
distribution of primordial discs, as would be expected in any 
model in which disc holes develop at a later evolutionary 
stage. However, we find a significant minority of transition 
discs at all mm fluxes including very high mm fluxes similar 
to the incidence of transition discs in the u pper quartile of 
the pr imordial disc mm distribution noted bv lAndrews et all 
|201ll ). 

We explore the hypothesis that the mm bright transi- 
tion discs might have a distinct physical origin by analysing 
the properties of two sub-samples of transition discs, split at 
the median of the primordial disc mm distribution. Virtually 
all the mm bright transition discs have high accretion rates 
and large inner holes. The mm faint sub-sample has a range 
of accretion rates and systematically lower hole sizes and 
broadly occupies the region of the accretion rate, hole size 
plane that is consistent with X-ray photoevaporation mod- 
els. The two sub-samples differ in accretion rate and inner 
hole size distributions at the > 4a level. Furthermore, early 
type stars are over represented in the mm bright sample 
compared to the primordial and faint mm transition discs. 
Nevertheless, the difference in properties cannot be solely 
ascribed to differences in host stars. A flat tail of mm bright 
sources are found at all spectral types and these have sys- 
tematically larger holes than mm faint objects of the same 
spectral type. Within the mm bright sub-sample there is a 
strong correlation between inner hole size and spectral type 
(uncertainties in hole sizes for many of the mm faint objects 
do not allow us to assess this correlation within the mm faint 
sub-sample). 

Finally, we ask whether there are any obvious selec- 
tion effects that drive the observed association between mm 
bright sources and large holes. We believe there are no fac- 
tors that would militate against the identification of small 
holes in discs that are mm bright. Moreover, since the evi- 
dence of even the largest inner holes in the sample is based 
on spectral dips at < 30(im, then the detection of such a 
hole would not be undermined by a low flux at 1.3mm. We 
therefore do not believe that it is a detectability issue that 
is driving this association. We thus need to ask a) why mm 
bright discs apparently do not form small holes in small 
and large dust grains and b) is there a population of discs 
with such large holes that they are classified as disc-less 
sources on the basis of their near-infrared to mm SEDs? 
This latter issue was raised by Owen et al (2011) although 
their subsequent work on Xray photoevaporation (Owen et 
al 2012) suggested that the process of 'thermal sweeping' 
should rapidly disperse such low surface density relic discs. 
We also notes that some discs may manifest holes in mm im- 
ages while still showing high levels high levels of NIR excess 
as hinted by Andrews et al. (2011), whether one would class 
these discs within the same existing 'transition' disc frame- 
work is unclear and certainly a matter for future debate. 

We stress that these findings are necessarily suggestive 
rather than conclusive and underline avenues of future inves- 
tigation. Two obvious shortcomings are i) the lack of fitted 
inner hole radii in a number of transition discs in the liter- 
ature and ii) the lack of good control samples for the mm 
properties of primordial discs at a range of spectral types. 
Clearly ALMA will assist in the latter respect. 
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